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We discuss several novel types of multi-component (temporal and spatial) envelope solitary waves 
that appear in fiber and waveguide nonlinear optics. In particular, we describe multi-channel solitary 
waves in bit-parallel- wavelength fiber transmission systems for high performance computer networks, 
multi-colour parametric spatial solitary waves due to cascaded nonlinearities of quadratic materials, 
and quasiperiodic envelope solitons due to quasi-phase-matching in Fibonacci optical superlattices. 



I. INTRODUCTION 

Rapid progress in the design and manufacture of op- 
tical fiber systems is a result of worldwide demand for 
ultra-high bit-rate optical communications. This ex- 
plains the growing interest of the soliton community 
in soliton-based optical fiber communication systems. 
This area of research was considerably advanced in re- 
cent years The most remarkable results include 
the application of the concept of the dispersion man- 
agement to temporal optical solitons and soliton-based 
optical transmission systems, and the discovery of the 
so-called dispersion managed soliton. High-speed opti- 
cal communications require effective components such 
as high-performance broadband computer networks that 
can be developed by employing the concept of the bit- 
parallel-wavelength (BPW) pulse transmission that of- 
fers many of the advantages of both parallel fiber ribbon 
cable and conventional wavelength-division-multiplexing 
(WDM) systems §. 

Expanding development in the study of the soliton 
fiber systems has been observed in parallel with impres- 
sive research on their spatial counterparts, optical self- 
trapped beams or spatial optical solitons. One of the 
key concepts in this field came from the theory of multi- 
frequency wave mixing and cascaded nonlinearities where 
a nonlinear phase shift is produced as a result of the para- 
metric wave interaction In all such systems, the 
nonlinear interaction between the waves of two (or more) 
frequencies is the major physical effect that can support 
coupled-mode multi-frequency solitary waves. 

The examples of temporal and spatial solitons men- 
tioned above have one common feature: they involve 
the study of solitary waves in multi-component nonlinear 
models. The main purpose of this paper is to overview 
several different physical examples of multi-mode and/or 
multi-frequency solitary waves that occur for the pulse or 
beam propagation in nonlinear optical fibers and waveg- 
uides. For these purposes, we select three different 
cases: multi-wavelength solitary waves in bit-parallel- 
wavelength optical fiber links, multi-colour spatial soli- 
tons due to multistep cascading in optical waveguides 
with quadratic nonlinearities, and quasiperiodic solitons 
in the Fibonacci superlattices. We believe these examples 



display both the diversity and richness of the multi-mode 
soliton systems, and they will allow further progress to be 
made in the study of nonlinear waves in multi-component 
nonintegrable physical models. 



II. TEMPORAL AND SPATIAL SOLITONS 

Because the phenomenon of the long-distance prop- 
agation of temporal optical solitons in optical fibers 0| 
is known to a much broader community of researchers 
in optics and nonlinear physics, first we emphasize the 
difference between temporal and spatial solitary waves. 
Indeed, for a long time stationary beam propagation in 
planar waveguides has been considered somewhat simi- 
lar to the pulse propagation in fibers. This approach is 
based on the so-called spatio-temporal analogy in wave 
propagation, meaning that the propagation coordinate z 
is treated as the evolution variable and the spatial beam 
profile along the transverse direction in waveguides, is 
similar to the temporal pulse profile in fibers. This anal- 
ogy is based on a simple notion that both beam evolution 
and pulse propagation can be described by the cubic non- 
linear Schrodinger (NLS) equation. 

However, contrary to the widely accepted opinion, 
there is a crucial difference between temporal and spa- 
tial solitons. Indeed, in the case of the nonstationary 
pulse propagation in fibers, the operation wavelength is 
usually selected near the zero point of the group- velocity 
dispersion. This means that the absolute value of the 
fiber dispersion is small enough to be compensated by 
a weak nonlinearity such as that produced by the (very 
weak) Kerr effect in optical fibers which leads to a rela- 
tive nonlinearity-induced change in the refractive index. 
Therefore, nonlinearity in such systems is always weak 
and it should be well modeled by a cubic NLS equa- 
tion which is known to be integrable by means of the 
inverse-scattering technique. However, for very short 
(e.g., fs) pulses the cubic NLS equation describing the 
long-distance propagation of pulses should be corrected 
to include additional terms that would account for such 
effects as higher-order dispersion, Raman scattering, etc. 
All such corrections can be taken into account with the 
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help of the perturbation theory Thus, in fibers non- 
linear effects are weak and they become important only 
when dispersion is small (near the zero-dispersion point) 
affecting the pulse propagation over large distances (of 
order of hundreds of meters or even kilometers). The sit- 
uation changes dramatically when we consider the prop- 
agation of multi-wavelength pulses with almost equal 
group velocities. The corresponding model is described 
by a nonintegrable and rather complicated system of cou- 
pled NLS equations, which we briefly discuss below. 

In contrary to the pulse propagation in optical fibers, 
the physics underlying the stationary beam propagation 
in planar waveguides and bulk media is different. In 
this case an optical beam is generated by a continuous 
wave (CW) source and it is time independent. However, 
when the beam evolves with the propagation distance z, 
it diffracts in the transverse spatial directions. Then, a 
nonlinear change in the refractive index should compen- 
sate for the beam spreading caused by diffraction which is 
not a small effect. That is why to observe spatial solitons 
as self-trapped optical beams, much larger nonlinearities 
are usually required, and very often such nonlinearities 
are not of the Kerr type (e.g. they saturate at higher 
intensities). This leads to the models of generalized non- 
linearities with the properties of solitary waves different 
from those described by the integrable cubic NLS equa- 
tion. Propagation distances involved in the phenomenon 
of the beam self-focusing and spatial soliton propagation 
are of the order of millimeters or centimeters. To achieve 
such large nonlinearities, one needs to use the optical 
materials with large nonlinearity-induced refractive in- 
dex. One of the possible way to overcome this difhculty 
is to use the so-called cascaded nonlinearities of noncen- 
trosymmetric optical materials where nonlinear effects 
are accumulated due to parametric wave interaction un- 
der the condition of the wave phase matching. Such para- 
metric wave-mixing effects generate novel classes of spa- 
tial optical solitons where resonant parametric coupling 
between the envelopes of two (or more) beams of differ- 
ent frequencies supports stable spatially localised waves 
even in a bulk medium (see details in Ref. |3|). It is this 
kind of multi-component solitary waves that we discuss 
below. 



III. BIT-PARALLEL- WAVELENGTH SOLITONS 

A growing demand for high-speed computer commu- 
nications requires an effective and inexpensive computer 
interconnection. One attractive alternative to the con- 
ventional WDM systems is BPW single-mode fiber op- 
tics links for very high bandwidth computer communi- 
cations They differ from the WDM schemes in that 
no parallel to serial conversion is necessary, and paral- 
lel pulses are launched simultaneously on different wave- 
lengths. 

When the pulses of different wavelengths are transmit- 



ted simultaneously, the cross-phase modulation can pro- 
duce an interesting pulse shepherding effect when a 
strong ("shepherd") pulse enables the manipulation and 
control of pulses co-propagating on different wavelengths 
in a multi-channel optical fiber link. 

To describe the simultaneous transmission of N differ- 
ent wavelengths in a nonlinear optical fiber, we follow the 
standard derivation Q and obtain a system of N coupled 
nonlinear Schrodinger (NLS) equations (0 < j < N — 1): 
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where, for the j— th wave, Aj^z^t) is the slowly varying 
envelope, Vgj and /32j" are the group velocity and group- 
velocity dispersion, respectively, and the nonlinear coef- 
ficients Xj characterize the Kerr effect. Equations (^ do 
not include the fiber loss, since the fiber lengths involved 
in bit-parallel links are only a small fraction of the atten- 
uation length. 

We measure the variables in the units of the central 
wavelength channel (say, j — Q), and obtain the follow- 
ing normalized system of the coupled NLS equations. 




(2) 



Uj — 0, 



is the incident optical power 



where Uj = 

in the central channel, aj = (/32j/|/32o|), Ij = Xj/xo, so 
that Q!o = 7o = 1- For the operating wavelengths spaced 
1 nm apart within the band 1530 1560 nm, the coeffi- 
cients Uj and 7j are different but close to 1. Initially, in 
Eq. (^, we omit the mode walk-off effect described by 
the parameters 5j = v~q — (so that 6q — 0). This 
effect will be analysed later in this section. 

To analyze the nonlinear modes, i.e. localized states 
of the BPW model (||), we look for stationary solutions 
in the form, 



(3) 



and therefore obtain the system of equations for the nor- 
malized mode amplitudes. 



1 d^UQ 
2~dt^ 



N-l \ ^ 

|W0P + 2 ^ \Un\^ I Uq = -Uq, 
n=l / 
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2 " dt^ 



- In I + 2 ^ j M„ = XnUn, 



where n=l,2,...,A — 1, the amplitudes and time are 
measured in the units of y/2f3o and {2(3q)~^^^, respec- 
tively, and A„ = /3n/2/3o. 
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System has exact analytical solutions for N cou- 
pled components, the so-called BPW solitons. Indeed, 
looking for solutions in the form uo{t) — Uq sechi. 
Unit) = Un sech i, we obtain the constraint A„ = Q!„/2, 
and a system of N coupled algebraic equations for the 
wave amplitudes, 



n=l 



1, u' 



In a special symmetric case, we take a„ = 7„ = 1, and 
the solution of those equations is simple [Q: 
Uo = Un = U, = [l + 2(7V-l)]-i/2. 

Analytical solutions can also be obtained in the lin- 
ear limits when the central frequency pulse (at ti = 0) is 
large. Then, linearizing Eqs. (^) for small |u„| <^ |uo|, we 
obtain a decoupled nonlinear equation for mq and TV — 1 
decoupled linear equations for u„. Each of the latter 
possess a localized solution provided A„ — \n \ where 
Ak°^ = (a„/8)[l - v^l -t- 16(7„/a„)]^ In this limit the 
central soliton pulse mq ("shepherd pulse") can be con- 
sidered as inducing an effective waveguide that supports 
a fundamental mode w„ with the corresponding cutoff 
aI^'' . Since, by definition, the parameters q;„ and 7„ are 
close to 1, we can verify that the soliton-induced waveg- 
uide supports maximum of two modes (fundamental and 
the first excited one). This is an important physical re- 
sult that explains the effective robustness of the pulse 
guidance by the shepherding pulse. 
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FIG. 1. Bifurcation diagram (top) and two examples of 
the stationary BPW soliton solutions for A*' = 4. 

To demonstrate a number of unique properties of the 
multi-channel BPW solitons, we consider the case = 4 
in more details. A comprehensive discussion of the case 
N = 2 can be found in the preprint jsj . We select the fol- 
lowing set of the normalized parameters: ao.s = 7o,3 = 1, 



Q!i = 7i = 1.1, and 0.2 = ^2 = 1.05. Solitary waves of 
this four-wavelength BPW system can be found numeri- 
cally as localized solutions of Eqs. (^). Figure |] presents 
the lowest-order families of such localized solutions. In 
general, they are characterized by A — 1 parameters, but 
we can capture the characteristic features by presenting 
power dependencies along the line A = Ai = A2 = . • . 
in the parameter space {A„}. The power of the central- 
wavelength component (n = 0) does not depend on A 
(straight line Pq = 2). Thin dashed, dotted, and dash- 
dotted curves correspond to the three separate single- 
mode solitons of the multi-channel BPW system, (1), 
(2), and (3), respectively, shown with the corresponding 
branches of (O-f-l), (0+2), and (0-f3) two-mode solitons. 
The latter curves start off from the bifurcation points on 
the Mo branch at A^°^ , Aj"^ , and Ag"^ , respectively. Close 
separation of those curves is the result of closeness of the 
parameters a„ and 7„ for n = 1,2,3. 

Thick solid curves in Fig. ^ correspond to the two- 
(1-1-2) and three-mode (O-t-l-l-2) localized solutions. 
The latter solutions bifurcate and give birth to four- 
wavelength solitons (0-I-1-I-2-I-3) (branch A-B). Two ex- 
amples of such four-wave composite solitons are shown 
in Fig. |l| (bottom row). The solution B is close to an 
exact sech-type solution at A = 0.5 (described above) for 
A^ = 4, whereas the solution A is close to that approx- 
imately described in the linear limit in the vicinity of a 
bifurcation point. Importantly, for different values of the 
parameters (a„,7„), the uppermost bifurcation point for 
this branch (open circle in Fig. |^) is not predicted by 
a simple linear theory and, due to the nonlinear mode 
coupling, it gets shifted from the branch of the central- 
wavelength soliton (straight line) to a two-mode branch 
(0+1-^2) (thick solid curve). 

As a result, if we start on the right end of the horizontal 
branch and follow the lowest branches of the total power 
P(A) in Fig. |l|, we pass the following sequence of the 
soliton families and bifurcation points: (0) ^ (0 + 1) ^ 
(0+1+2) ^ (0+1+2+3) ^ (1+2+3) (2+3) -> (3). If 
the modal parameters are selected closer to each other, 
the first two links of the bifurcation cascade disappear 
(i.e. the corresponding bifurcation points merge), and 
the four-mode soliton bifurcates directly from the central- 
wavelength pulse, as predicted by the linear theory. Note 
however that the sequence and location of the bifurcation 
points is a function of the cross-section of the parameter 
space {A„}, and the results presented above correspond 
to the choice A = Ai = A2 = . . .. 

The qualitative picture of the cascading bifurcations 
preserves for other values of A^. In particular, near the bi- 
furcation point a mixed-mode soliton corresponds to the 
localized modes guided by the central-wavelength soliton 
(shepherd) pulse. The existence of such soliton solutions 
is a key concept of BPW transmission when the data are 
launched in parallel carrying a desirable set of bits of in- 
formation, all guided by the shepherd pulse at a selected 
wavelength. 



3 



Effects of the walk-off on the muhi-channel BPW soh- 
tons seems to be most dangerous for the pulse alignment 
in the parallel links. For nearly equal soliton components, 
it was shown long time ago that nonlinearity can pro- 
vide an effective trapping mechanism to keep the pulses 
together. For the shepherding effect, the correspond- 
ing numerical simulations are presented in Figs, ^(a-d) 
for the four-channel BPW system. Initially, we launch 
a composite four-mode soliton as an unperturbed solu- 
tion A [see Fig. |^ of Eqs. (||), without walk-off and cen- 
tered a,t t — 0. When this solution evolves along the 
propagation direction z in the presence of small to mod- 
erate relative walk-off (5„ ^ for n ^ 0), its compo- 
nents remain strongly localized and mutually trapped 
[Fig. ||(a,b)], whereas it loses some energy into radia- 
tion for much larger values of the relative mode walk-off 
[Fig. |(c,d)]. 



lowed by reconstruction of the fundamental wave through 
the down-conversion frequency mixing (DFM) process 
{2uj — uj ~ oj). These two processes are governed by one 
phase-matched interaction and they differ only in the di- 
rection of power conversion. 

The idea to explore more than one simultaneous 
nearly phase-matched process, or double-phase-matched 
(DPM) wave interaction^ became attractive only re- 
cently jl^Jll[|, for the purposes of all-optical transistors, 
enhanced nonlinearity-induced phase shifts, and polar- 
ization switching. In particular, it was shown pT| that 
multistep cascading can be achieved by two second-order 
nonlinear cascading processes, SHG and sum-frequency 
mixing (SFM), and these two processes can also support a 
novel class of multi-colour parametric solitons [|l2| , briefly 
discussed below. 
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FIG. 2. Influence of the mode walk-ofT on the shep- 
herding effect for the input pulse A of Fig. |l|. Shown 
are: output pulse at the propagation distance z = 50 
for (a) 5i = 0.45,(52 = -0.35,53 = 0.25, and 
(c) (5i = J2 = ^3 = 0.9; (b,d) — evolution of the position 
of maximum for all pulse constituents. Solid line — maxima 
of the shepherd pulse. 



IV. PARAMETRIC OPTICAL SOLITONS DUE 
TO MULTISTEP CASCADING 

A. Concept of multistep cascading 

Recent progress in the study of cascading effects in 
optical materials with quadratic (second-order or x^^-*) 
nonlinear response has offered new opportunities for all- 
optical processing, optical communications, and optical 
solitons 1^,^. Most of the studies of cascading effects 
employ parametric wave mixing processes with a sin- 
gle phase- matching and, as a result, two-step cascad- 
ing iQj. For example, the two-step cascading associated 
with type I second-harmonic generation (SHG) includes 
the generation of the second harmonic [uj + u) = 2uj) fol- 



B. Multistep cascading solitons 

To introduce the simplest model of multistep cascad- 
ing, we consider the fundamental beam with frequency uj 
entering a noncentrosymmetric nonlinear medium with a 
quadratic response. As a first step, the second-harmonic 
wave with frequency 2uj is generated via the SHG pro- 
cess. As a second step, we expect the generation of 
higher order harmonics due to SFM, for example, a 
third harmonic (uj -t- 2a; = 3w) or even fourth harmonic 
{2uj-\-2lli = Alo) When both such processes are nearly 
phase matched, they can lead, via down-conversion, to a 
large nonlinear phase shift of the fundamental wave [^ . 
Additionally, the multistep cascading can support a novel 
type 0/ three-wave spatial solitary waves in a diffractive 
X^"^^ nonlinear medium, multistep cascading solitons. 

We start our analysis with the reduced amplitude equa- 
tions derived in the slowly varying envelope approxima- 
tion with the assumption of zero absorption of all inter- 
acting waves (see, e.g., Ref. (Q). Introducing the effect 
of diffraction in a slab waveguide geometry, we obtain 
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X-^A^Axe 



0, 



where xi,2 = 2k\a\^2, X3 = 6fcio-3, and X4,5 = 4A;i(T4,5, 
and the nonlinear coupling coefficients crj, are propor- 
tional to the elements of the second-order susceptibility 
tensor which we assume to satisfy the following relations 
(no dispersion), 0-3 = Sci, (T2 = and cr4 = 2a \. 

In Eqs. (||), A\,Ai and Az are the complex electric 
field envelopes of the fundamental harmonic (FH), second 
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harmonic (SH), and third harmonic (TH), respectively, 
Afc2 = 2ki — k2 is the wavevector mismatch for the SHG 
process, and Ak^ — ki + k2 ~ k^ is the wavevector mis- 
match for the SFM process. The subscripts '1' denote 
the FH wave, the subscripts '2' denote the SH wave, and 
the subscripts '3', the TH wave. FoUowing the technique 
earher employed in Refs. we look for stationary so- 
lutions of Eq. (|^) and introduce the normalised envelope 
w{z,x), v{z,x), and u{z,x) according to the relations. 
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where Ak = Ak2 + Ak^ 
z z/ 13 and x xj \J2fik\ 
of coupled equations. 



Renormalising the variables as 
we finally obtain a system 



dw 

oz 
dv 
dz 
du 
dz 



dx^ 
dx"^ 
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— w + w*v + v*u = 0, 

— av ~V —w^ + w*u = 0, 

— aiu + xvw = 0, 



(7) 



where a = 2(2/3 + Afe)//^ and ai = 3(3/3 -I- Afc)//3 
are two dimensionless parameters that characterise the 
nonlinear phase matching between the parametrically 
interacting waves. Dimensionless material parameter 
X = X1X3/X2 = 9(cri/CT2)^ depends on the type of 
phase matching, and it can take different values of or- 
der of one. For example, when both SHG and SFM 
are due to quasi-phase matching (QPM), we have (Tj = 
(2/7rTO)(7r/Aini)x(^)[w;(4- jV;-(3- j>], where j = 
1,2. Then, for the first-order (m — 1) QPM processes 
(see, e.g., Ref. we have ai = 0-2, and therefore 

X = 9. When SFM is due to the third-order QPM pro- 
cess (see, e.g., Ref. we should take cti = (72/3, and 
therefore x = 1. At last, when SFM is the fifth-order 
QPM process, we have cti — <72/5 and x — 9/25. 

Dimensionless equations (^ present a fundamental 
model for three-wave multistep cascading solitons in 
the absence of walk-off. Additionally to the type I 
SHG solitons (see, e.g., Refs Q), the multistep cas- 
cading solitons involve the phase-matched SFM inter- 
action (uj + 2uj = 3a;) that generates a third harmonic 
wave. Two-parameter family of localised solutions con- 
sists of three mutually coupled waves. It is interesting to 
note that, similar to the case of nondegenerate three- 
wave mixing iQ, Eqs. (|^) possess an exact solution. 
To find it, we make a substitution w = sech^ (r]x) , 
V = Wo sech^(?7a:;) and u — uq sech^ (rjx) , and obtain un- 
known parameters from the following algebraic equations 



3 + 'ix^o ' 



4x^0 + 6wo = 9, Mo = xX'f^owo, (8) 



valid for V — ^ ^^^d a — ai — 1. Equations (g) have 
two solutions corresponding to positive and negative val- 
ues of the amplitude (wq). This indicates a possibility 
of multi-valued solutions, even within the class of exact 
solutions. 
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FIG. 3. Examples of three-wave solitary waves of the 
model (^) for (a) a — 0.05, Qi — 5, and (b) a = 5, ai = 0.35. 

In general, three- wave solitons of Eqs. (Q) can be found 
only numerically. Figures |^(a) and |^(b) present two ex- 
amples of solitary waves for different sets of the mismatch 
parameters a and ai. When ai ^ I [see Fig. ||(a)], 
which corresponds to an unmatched SFM process, the 
amplitude of the third harmonic is small, and it vanishes 
for Qfi — > 00. 
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a, 

FIG. 4. Two branches of multistep cascading solitons 
shown as the total soliton power P vs. the parameter ai 
for Q = 1 and x ~ ^- Filled circles show the analytical so- 
lutions. The lower branch approaches a family of two-wave 
quadratic solitons (for ai — » 00) shown by a dashed line. 

To summarise different types of three-wave solitary 
waves, in Fig. || we plot the dependence of the total soli- 
ton power defined as 



P^J^ dx (^\w\^ + A\v\^ + ^\u\^'^ 



(9) 



on the mismatch parameter ai, for fixed a = 1. It is 
clearly seen that for some values of ai (including the ex- 
act solution at ai = 1 shown by two filled circles), there 
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exist two different branches of three-wave solitary waves, 
and only one of those branches approaches, for large val- 
ues of ai, a family of two-wave solitons of the cascad- 
ing limit (Fig. ^, dashed). The slope of the branches 
changes from negative (for small ai) to positive (for large 
ai), indicating a possible change of the soliton stability. 
However, the detailed analysis of the soliton stability is 
beyond the scope of this paper (see, e.g., Refs. [p^ , p^ ). 



C. Parametric soliton-induced waveguides 

Another type of multistep cascading parametric pro- 
cesses which involve only two frequencies, i.e. two-colour 
multistep cascading, can occur due to the vectorial in- 
teraction of waves with different polarization. We de- 
note two orthogonal polarization components of the fun- 
damental harmonic (FH) wave {u>i = w) as A and B, 
and two orthogonal polarizations of the second harmonic 
(SH) wave {uj2 = 2a;), as S and T. Then, a simple mul- 
tistep cascading process consists of the following steps. 
First, the FH wave A generates the SH wave S via type 
I SHG process. Then, by down-conversion SA-B, the or- 
thogonal FH wave B is generated. At last, the initial 
FH wave A is reconstructed by the processes SB-A or 
AB-S, SA-A. Two principal second-order processes AA- 
S and AB-S correspond to two different components of 
the x^^^ susceptibility tensor, thus introducing additional 
degrees of freedom into the parametric interaction. Dif- 
ferent cases of such type of multistep cascading processes 
are summarized in Table ^ 

TABLE I. Two-frequency multistep cascading processes 





Principal 


Equivalent 


(a) 


(AA-S, AB-S) 


(BB-S, AB-S); (AA-T, AB-T) 
(BB-T, AB-T) 


(b) 


(AA-S, AB-T) 


(BB-S, AB-T); (AA-T, AB-S) 
(BB-T, AB-S) 


(c) 


(AA-S, BB-S) 


(AA-T, BB-T) 


(d) 


(AA-S, AA-T) 


(BB-S, BB-T) 



To demonstrate some of the unique properties of 
the multistep cascading, we discuss here how it can 
be employed for soliton-induced waveguiding effects in 
quadratic media. For this purpose, we consider a model 
of two-frequency multistep cascading described by the 
principal DPM process (c) (see Table | above) in the pla- 
nar slab-waveguide geometry. Using the slowly varying 
envelope approximation with the assumption of zero ab- 
sorption of all interacting waves, we obtain 
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(10) 



where X\-,i = 2fci(Ti_2, the nonlinear coupling coefficients 
CTfc are proportional to the elements of the second-order 
susceptibility tensor, and Afci and Afc2 are the corre- 
sponding wave- vector mismatch parameters. 

To simplify the system (p^, we look for its station- 
ary solutions and introduce the normalized envelopes 
M, u, and w according to the following relations, A = 
7iM exp(i/3z — ^A/ciz), B = 72^ exp(i/?z — ^l\kiz\ 
and S — 73W exp(2i/?z), where 7^""'^ — 2xiXq, — 
2x§(xiX2)^^^, and 7;^^ = Xi^o: and the longitudinal 
and transverse coordinates are measured in the units of 
zo — {0 — and xq = (zo/2fci)^/^, respectively. 

Then, we obtain a system of normalized equations. 



■du d^u , 

+ "5^ -u + u w = 0, 
oz ox'' 

.dv d^v 



2i 



dz dx^ 
dw d^w 



dz 



dx^ 



aiv + XV* w — 0, 
— aw H — (u^ + V 



(11) 



0, 



where x = (X2/X1), ai = (/? - Afc2/2)(/3 - AA:i/2)-i, 
and a = 4/? {f3-Aki/2)-\ 

First of all, we notice that for v = (or, similarly, 
u = 0), the dimensionless Eqs. ([ll| ) reduce to the corre- 
sponding model for the two-step cascading due to type I 
SHG discussed earlier ||,^, and its stationary solutions 
are defined by the equations for real u and w, 



d u 
dx^ 



u + uw = 0, 
1 



(12) 



aw ■ 



-u 



0, 



that possess a one-parameter family of two-wave local- 
ized solutions (mq , ) found earlier numerically for any 
Q ^ 1, and also known analytically for a = 1, Uf){x) = 
(3/%/2)sech2(a;/2) = V2wq{x) (see Ref. @). 

Then, in the small-amplitude approximation, the equa- 
tion for real orthogonally polarized FH wave v can be 
treated as an eigenvalue problem for an effective waveg- 
uide created by the SH field wo{x), 



dx'^ 



[xwaix) - ai]v = 0. 



(13) 



Therefore, an additional parametric process allows to 
propagate a probe beam of one polarization in an effec- 
tive waveguide created by a two-wave spatial soliton in 
a quadratic medium with FH component of another po- 
larization. However, this type of waveguide is different 
from what has been studied for Kerr-like solitons because 
it is coupled parametrically to the guided modes and, as a 
result, the physical picture of the guided modes is valid, 
rigorously speaking, only in the case of stationary phase- 
matched beams. As a result, the stability of the corre- 
sponding waveguide and localized modes of the orthog- 
onal polarization it guides is a key issue. In particular. 
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the waveguide itself (i.e. two-wave parametric soliton) 
becomes unstable for a < acr ~ 0.2 |lq] . 

In order to find the guided modes of the parametric 
waveguide created by a two- wave quadratic soliton, we 
have to solve Eq. (O) where the exact solution ^0(2^) is 
to be found numerically. Then, to address this problem 
analytically, approximate solutions can be used, such as 
those found with the help of the variational method . 
However, the different types of the variational ansatz 
used do not provide a very good approximation for the 
soliton profile at all a. For our eigenvalue problem 
( p^ ) , the function wq [x) defines parameters of the guided 
modes and, in order to obtain accurate results, it should 
be calculated as close as possible to the exact solutions 
found numerically. To resolve this difficulty, below we 
suggest a novel "almost exact" solution that would al- 
low to solve analytically many of the problems involving 
quadratic solitons, including the eigenvalue problem (13). 

First, we notice that from the exact solution at a = 1 
and the asymptotic result for large a, w ~ u^/{2a) « 
2 sech^(a;), it follows that the SH component Wo{x) of 
Eqs. ( p^ remains almost self-similar for a > 1. Thus, we 
look for the SH field in the form wo{x) = Wm sech^(a;/p), 
where Wm and p are to be defined. The solution for uq(x) 
should be consistent with this choice of the shape for SH, 
and it is defined by the first (linear for u) equation of the 
system (^2|). Therefore, we can take u in the form of the 
lowest guided mode, uq{x) = sech^(a;/p), that corre- 
sponds to an effective waveguide {x) . By matching the 
asymptotics of these trial functions with those defined di- 
rectly from Eqs. (^2|) at small and large x, we obtain the 
following solution, 

uq{x) = UmSech^(x/p), wo{x) = WmSech^ {x / p) , (14) 



1 



4(w™ - 1)' 



(t«m-l) (Wm-1) (2-W„) 

Here, the third relation allows us to find Wm for arbi- 
trary a as a solution of a cubic equation, and then to 
find all other parameters as functions of a. For mis- 
matches in the interval < a < -|-cxd, the parameter val- 
ues change monotonically in the regions: < u„i < +00, 
1 < Wm < 2, and -\-oo > p > 1. It is really amazing that 
the analytical solution ( p^ ) , ( p^ provides an excellent ap- 
proximation for the profiles of the two-wave parametric 
solitons found numerically, with the relative errors not 
exceeding l%-3% for stable solitons (e.g. when a > aci). 
As a matter of fact, we can treat Eqs. (|lj) and ([l^) as an 
approximate scaling transformation of the family of two- 
wave bright solitons. Moreover, this solution allows us to 
capture some remarkable internal similarities and distinc- 
tions between the solitons existing in different types of 
nonlinear media. In particular, as follows from Eqs. (|lj) 
and (|l5|), the FH component and the self-consistent ef- 
fective waveguide (created by the SH field) have approxi- 
mately the same stationary transverse profiles as for one- 
component solitons in a Kerr-like medium with power- 
law nonlinear response pll]. For a = 1 (p = 2) and a ^ 1 



(p = 1) our general expressions reduce to the known an- 
alytical solutions, and the FH profile is exactly the same 
as that for solitons in quadratic and cubic Kerr media, 
respectively. On the other hand, the strength of self- 
action for quadratic solitons depends on the normalized 
phase mismatch a and, in general, the beam dynamics 
for parametric wave mixing can be very different from 
that observed in Kerr-type media. 

Now, the eigenvalue problem ( |l3|) can be readily solved 
analytically. The eigenmode cutoff values are defined 
by the parameter ai that takes one of the discrete val- 
ues, a^"^ = (s - where s = -(1/2) + [(1/4) -|- 
WmXP^]^/^. Number n stands for the mode order (n = 
0, 1, . . .), and the localized solutions are possible provided 
n < s. The profiles of the corresponding guided modes 
are 

v„{x) = Fsech''""(a;/p)iJ(-n,2s-n + l,s-n+l;C/2), 

where C = 1 — tanh(a;/p), H is the hypergeometric func- 
tion, and V is the mode's amplitude which cannot be 
determined within the framework of the linear analysis. 




(b) 


n = / 




1 



FIG. 5. Cutoff eigenvalues ol\^^ of the guided modes shown 
as (a) functions of a at x = 2, and (b) functions of x at a = 4. 
Dashed lines correspond to the intersection of the plots in the 
parameter space (ol,\). 

According to these results, a two- wave parametric soli- 
ton creates, a multi-mode waveguide and larger number 
of the guided modes can be observed for smaller a. Fig- 
ures ||(a,b) show the dependence of the mode cutoff val- 



ues cx'^\a) for a fixed Xi and oi\\x) for a fixed a, re- 



spectively. 



For the case x 



1, the dependence has a 
simple form: aY'''{a) = [1 — n{wm — 1)]^- 

Because a two-wave soliton creates an induced waveg- 
uide parametrically coupled to its guided modes of the 
orthogonal polarization, the dynamics of the guided 
modes may differ drastically from that of conventional 
waveguides based on the Kerr-type nonlinearities. Fig- 
ures |6|(a-d) show two examples of the evolution of guided 
modes. In the first example [see Fig. ^(a-c)], a weak fun- 
damental mode is amplified via parametric interaction 
with a soliton waveguide, and the mode experiences a 
strong power exchange with the orthogonally polarized 
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FH component through the SH field. This process is ac- 
companied by only a weak deformation of the induced 
waveguide [see Fig. |(a) - dotted curve!. The resulting 
effect can be interpreted as a power exchange between 
two guided modes of orthogonal polarizations in a waveg- 
uide created by the SH field. In the second example, the 
propagation is stable [see Fig. ||(d)]. 

When all the fields in Eq. ( pi] ) are not small, i.e. the 
small-amplitude approximation is no longer valid, the 
profiles of the three-component solitons should be found 
numerically. However, some of the lowest-order states 
can be calculated approximately using the approach of 
the "almost exact" solution (|lj),(|l^) described above, 
which is presented in detail elsewhere [^ , ^ . More- 
over, a number of the solutions and their families can 
be obtained in an explicit analytical form. For exam- 
ple, for ai = 1/4, there exist two families of three- 
component solitary waves for any a > 1, that de- 
scribe soliton branches starting at the bifurcation points 
ai — a^^^ at a = 1: (i) the soliton with a zero-order 
guided mode for x — 1/3: u{x) = (3/-\/2) sech^ (a;/2), 
v{x) = C2 sech (a;/2), w{x) = (3/2) sech^ (x/2), and 
(ii) the soliton with a first-order guided mode for x = 1: 
u{x) = ci sech^ (a:/2), v{x) = C2 sech^ (a::/2) sinh (a;/2), 
'w{x) = (3/2) sech^ (2^/2), where C2 — \/3 (a — 1) and 
ci = v'{9/2) + c^. 



For a practical realization of the DPM processes and 
the soliton waveguiding effects described above, we can 
suggest two general methods. The first method is based 
on the use of two commensurable periods of the quasi- 
phase-matched (QPM) periodic grating. Indeed, to 
achieve DPM, we can employ the first-order QPM for 
one parametric process, and the third-order QPM, for 
the other parametric process. Taking, as an example, 
the parameters for LiNbOs and AA-S {xx — z) and BB-S 
{zz — z) processes [^, we find two points for DPM at 
about 0.89 /im and 1.25 /^m. This means that a single 
QPM grating can provide simultaneous phase-matching 
for two parametric processes. For such a configuration, 
we obtain x ~ 1-92 or, interchanging the polarization 
components, x ~ 0.52. The second method to achieve 
the conditions of DPM processes is based on the idea 
of quasi-periodic QPM grating |2^j2^ . Specifically, Fi- 
bonacci optical superlattices provide an effective way to 
achieve phase-matching at several incommensurable pe- 
riods allowing multi-frequency harmonic generation in a 
single structure. We describe the properties of such struc- 
tures in the next section. 



V. ENVELOPE SOLITONS IN FIBONACCI 
SUPERLATTICES 

A. Incoherence and solitary waves 




6 12 18 2430 




FIG. 6. (a) Change of tlie normalized power in FH (it, 
solid) and SH (w, dotted) components, which initially consti- 
tute a two-wave soliton, and in the guided mode (w, dashed) 
at X = 2, demonstrating amplification of a guided wave. Evo- 
lution of the guided wave and effective waveguide (SH) is 
presented in plots (b) and (c), respectively, (d) Stationary 
propagation of a stable fundamental mode (x = 1). For all 
the plots a = 4, the initial amplitude is vo = 0.1, and qi 
corresponds to the bifurcation point. 



For many years, solitary waves have been considered 
as coherent localized modes of nonlinear systems, with 
particle-like dynamics quite dissimilar to the irregular 
and stochastic behavior observed for chaotic systems [^ . 
However, about 20 years ago Akira Hasegawa, while de- 
veloping a statistical description of the dynamics of an 
ensemble of plane waves in nonlinear strongly dispersive 
plasmas, suggested the concept of a localized envelope 
of random phase waves |2|]. Because of the relatively 
high powers required for generating self-localized random 
waves, this notion remained a theoretical curiosity until 
recently, when the possibility to generate spatial optical 
solitons by a partially incoherent source was discovered 
in a photorefractive medium [p9| . 

The concept of incoherent solitons can be compared 
with a different problem: the propagation of a soliton 
through a spatially disordered medium. Indeed, due to 
random scattering on defects, the phases of the individual 
components forming a soliton experience random fluctu- 
ations, and the soliton itself becomes partially incoher- 
ent in space and time. For a low-amplitude wave (linear 
regime) spatial incoherence is known to lead to a fast 
decay. As a result, the transmission coefficient vanishes 
exponentially with the length of the system, the phe- 
nomenon known as Anderson localization [ pO[ . However, 
for large amplitudes (nonlinear regime), when the non- 
linearity length is much smaller than the Anderson local- 
ization length, a soliton can propagate almost unchanged 
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through a disordered medium as predicted theoretically 
in 1990 [0 and recently verified experimentally |^ . 

These two important physical concepts, spatial self- 
trapping of light generated by an incoherent source in 
a homogeneous medium, and suppression of Anderson 
localization for large-amplitude waves in spatially disor- 
dered media, both result from the effect of strong nonlin- 
earity. When the nonlinearity is sufficiently strong it acts 
as an effective phase-locking mechanism by producing a 
large frequency shift of the different random-phase com- 
ponents, and thereby introducing an effective order into 
an incoherent wave packet, thus enabling the formation 
of localized structures. In other words, both phenomena 
correspond to the limit when the ratio of the nonlinearity 
length to the characteristic length of (spatial or tempo- 
ral) fluctuations is small. In the opposite limit, when this 
ratio is large, the wave propagation is practically linear. 

Below we show that, at least for aperiodic inhomoge- 
neous structures, solitary waves can exist in the inter- 
mediate regime in the form of quasiperiodic nonlinear 
localized modes. As an example, we consider SHG and 
nonlinear beam propagation in Fibonacci optical super- 
lattices, and demonstrate numerically the possibility of 
spatial self-trapping of quasiperiodic waves whose enve- 
lope amplitude varies quasiperiodically, while still main- 
taining a stable, well-defined spatially localized structure, 
a quasiperiodic envelope soliton. 



B. Quasi-phase-matching in optical superlattices 



We consider the interaction of a fundamental wave 
with the frequency lo (FH) and its SH in a slab waveguide 
with quadratic (or x''^'') nonlinearity. Assuming the x*-^' 
susceptibility to be modulated and the nonlinearity to be 
of the same order as diffraction, we write the dynamical 
equations in the form 



A DA A TiA DA 



.du 



.dw 



2d^ 



dz 4 dx"^ 



d{z)u*ve-'^' = 0, 
I- d{z)u^e"^' = 0, 



(16) 



where u{x, z) and w(x, z) are the slowly varying en- 
velopes of the FH and SH, respectively. The parame- 
ter (3 — Ak\kt^\xQ is proportional to the phase mismatch 
Afc = 2ki^ — k2uj, k^ and k2uj being the wave numbers 
at the two frequencies. The transverse coordinate x is 
measured in units of the input beam width xq, and the 
propagation distance z in units of the diffraction length 
Id = XQ\k;^\. The spatial modulation of the x*-^-* suscep- 
tibility is described by the quasi-phase-matching (QPM) 
grating function d{z). In the context of SHG, the QPM 
technique is an effective way to achieve phase matching, 
and it has been studied intensively [p3[. 




20 -10 60 HO 100 
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FIG. 7. (a) Slab waveguide with quasiperiodic QPM su- 
perlattice structure composed of building blocks A and B. (b) 
Numerically calculated amplitude spectrum of d{z). 

Here we consider a QPM grating produced by a 
quasiperiodic nonlinear optical superlattice. Quasiperi- 
odic optical superlattices, one-dimensional analogs of 
quasicrystals Q, are usually designed to study the ef- 
fect of Anderson localization in the linear regime of light 
propagation. For example, Gellermann et al. measured 
the optical transmission properties of quasiperiodic di- 
electric multilayer stacks of Si02 and Ti02 thin films and 
observed a strong suppression of the transmission 
For QPM gratings, a nonlinear quasiperiodic superlat- 
tice of LiTaOs, in which two antiparallel ferro-electric do- 
mains are arranged in a Fibonacci sequence, was recently 
fabricated by Zhu et al. , who measured multi-colour 
SHG with energy conversion efficiencies of ~ 5% — 20%. 
This quasiperiodic optical superlattice in LiTaOa can also 
be used for efficient direct third harmonic generation [p6| . 

The quasiperiodic QPM gratings have two building 
blocks A and B of the length I a and Ib, respectively, 
which are ordered in a Fibonacci sequence [Fig. 0(a)]. 
Each block has a domain of length Iai—'^ (^Si=1) with 
d=-\-l (shaded) and a domain of length Ia2—1{^ + v) 
[Ib2—1{^ — TT])] with d=—l (white). In the case of x^^-* 
nonlinear QPM superlattices this corresponds to positive 
and negative ferro-electric domains, respectively. The 
specific details of this type of Fibonacci optical superlat- 
tices can be found elsewhere |^^. For our simulations 
presented below we have chosen r]= 2(t — l)/(l-t-r^)= 
0.34, where r= (1 -I- -\/5)/2 is the so-called golden ra- 
tio. This means that the ratio of length scales is also the 
golden ratio, Ia/Ib= t. Furthermore, we have chosen 
^==0.1. 

The grating function d{z), which varies between -1-1 
and —1 according to the Fibonacci sequence, can be ex- 
panded in a Fourier series 



d{z) = dm,ne 



iG„ 



27r(m + nr) 
D '■ 



(17) 



where D=tIa + Ib=0-5'2 for the chosen parameter values. 
Hence the spectrum is composed of sums and differences 
of the basic wavenumbers ki=2'k/D and K2=2'itt/D. 
These components fill the whole Fourier space densely, 
since ki and K2 are incommensurate. Figure 0(b) shows 
the numerically calculated Fourier spectrum Gm,n- The 
lowest-order "Fibonacci modes" are clearly the most in- 
tense. 
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C. Quasiperiodic optical solitons 

To analyze the beam propagation and SHG in a 
quasiperiodic QPM grating one could simply average 
Eqs. (p^. To lowest order this approach always yields 
a system of equations with constant mean-value coeffi- 
cients, which does not allow to describe oscillations of 
the beam amplitude and phase. However, here we wish 
to go beyond the averaged equations and consider the 
rapid large-amplitude variations of the envelope func- 
tions. This can be done analytically for periodic QPM 
gratings However, for the quasiperiodic gratings we 
have to resolve to numerical simulations. 

Thus we have solved Eqs. ( |l6|) numerically with a 
second-order split-step routine. At the input of the crys- 
tal we excite the fundamental beam (corresponding to 
unseeded SHG) with a Gaussian profile, 



M(a;,0) = ^„e-^'/i°, w{x,0) = 0. 



(18) 



We consider the quasiperiodic QPM grating with match- 
ing to the peak at G2.3, i.e., /3=G2,3=82.25. First, we 
study the small-amplitude limit when a weak FH is in- 
jected with a low amplitude. Figures ||(a,b) show an 
example of the evolution of FH and SH in this effectively 
linear regime. As is clearly seem from Fig. ||(b) the SH 
wave is excited, but both beams eventually diffract. 

When the amplitude of the input beam exceeds a cer- 
tain threshold, self-focusing and localization should be 
observed for both harmonics. Figures ||(c,d) show an ex- 
ample of the evolution of a strong input FH beam, and its 
corresponding SH. Again the SH is generated, but now 
the nonlinearity is so strong that it leads to self-focusing 
and mutual self-trapping of the two fields, resulting in 
a spatially localized two-component soliton, despite the 
continuous scattering of the quasiperiodic QPM grating. 
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FIG. 8. (a) Diffraction of a weak FH beam with amplitude 
j4u=0.25 for /3=82.25. (c) Excitation of a quasiperiodic soli- 
ton by a FH beam with amplitude Au~5 for /3=82.25. (b,d) 
Corresponding SH components. 



It is important to notice that the two-component lo- 
calized beam created due to the self-trapping effect is 
quasiperiodic by itself. As a matter of fact, after an ini- 
tial transient its amplitude oscillates in phase with the 
quasiperiodic QPM modulation d{z). This is illustrated 
in Fig. ^, where we show in more detail the peak in- 
tensities in the asymptotic regime of the evolution. The 
oscillations shown in Fig. || are in phase with the oscilla- 
tions of the QPM grating d{z), and we indeed found that 
their spectra are similar. 
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FIG. 9. Amplitude oscillations of the quasiperiodic soli- 
ton. (a),(b) Close-ups of the peak intensity |u(0, z)] of the 
FH (black) and \w{0,z)\^ of the SH (grey). The Fibonacci 
building blocks A and B are indicated in (b) with d=l in 
grey regions, and d=—l in white regions. 



Our numerical results show that the quasiperiodic en- 
velope solitons can be generated for a broad range of the 
phase-mismatch /3. The amplitude and width of the soli- 
tons depend on the effective mismatch, which is the sepa- 
ration between (3 and the nearest strong peak Gm,n in the 
Fibonacci QPM grating spectrum [see Fig. |(b)]. Thus, 
low-amplitude broad solitons are excited for /3-values in 
between peaks, whereas high-amplitude narrow solitons 
are excited when (3 is close to a strong peak, as shown in 
Fig. |(c,d). 
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FIG. 10. (a) Transmission of the FH, |u(0, L)/m(0, 0)^, 
and the SH, |w(0, I/)/m(0, 0)^ vs. input peak intensity 
\u{0,0)f of the FH. (b) Output beam width vs. |u(0,0)|^ 
L=100, /?=Gi,2=50.83 (solid), /3=G2,4=101.66 (dashed), and 
/3=G2,3=82.25 (dot-dashed). 

To analyse in more detail the transition between the 
linear (diffraction) and nonlinear (self-trapping) regimes, 
we have made a series of careful numerical simula- 
. In Fig. 



tions 



10 we show the transmission coeffi- 
cients and the beam widths at the output of the crystal 
versus the intensity of the FH input beam, for a vari- 
ety of /?- values. These dependencies clearly illustrate the 
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universality of the generation of localised modes for vary- 
ing strength of nonlinearity, i.e. a quasiperiodic soliton is 
generated only for sufficiently high amplitudes. This is of 
course a general phenomenon also observed in many non- 
linear isotropic media. However, here the self- trapping 
occurs for quasiperiodic waves, with the quasiperiodicity 
being preserved in the variation of the amplitude of both 
components of the soliton. 



VI. CONCLUSION 

We have overviewed several important physical exam- 
ples of the multi-component solitary waves which appear 
due to multi-mode and/or multi- frequency coupling in 
nonlinear optical fibers and waveguides. We have de- 
scribed several types of such multi-component solitary 
waves, including: (i) multi-wavelength solitary waves 
in multi-channel bit-parallel-wavelength fiber transmis- 
sion systems, (ii) multi-colour parametric spatial solitary 
waves due to multistcp cascading in quadratic materi- 
als, and (iii) quasiperiodic envelope solitons in Fibonacci 
optical super lattices. These examples reveal some gen- 
eral features and properties of multi-component solitary 
waves in nonintegrable nonlinear models, also serving as 
a stepping stone for approaching other problems of the 
multi-mode soliton coupling and interaction. 
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